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ABSTRACT 

We report results from an intensive multiwavelength campaign on the 
intermediate-frequency-peaked BL Lacertae object W Com (z=0.102) during a 
strong outburst of very high energy gamma-ray emission in June 2008. The very 
high energy gamma-ray signal was detected by VERITAS on 2008 June 7-8 with 
a flux F(> 200 GeV) = (5.7 ± 0.6) x 10~ n cm _2 s _1 , about three times brighter 
than during the discovery of gamma-ray emission from W Com by VERITAS in 
2008 March. The initial detection of this flare by VERITAS at energies above 
200 GeV was followed by observations in high energy gamma-rays (AGILE, E 1 > 
100 MeV), and X-rays (Swift and XMM- Newton), and at UV, and ground-based 
optical and radio monitoring through the GASP-WEBT consortium and other 
observatories. Here we describe the multiwavelength data and derive the spec- 
tral energy distribution (SED) of the source from contemporaneous data taken 
throughout the flare. 

Subject headings: BL Lacertae objects: individual (W Com) - gamma rays: ob- 
servations 
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1. Introduction 

W Com (ON 231; z = 0.102) is a gamma-ray blazar classified as an intermediate- 



frequency-peaked BL Lac (IBL) object (INieppola et all 120061 ). based on the locations of its 



low-energy synchrotron peak and high-energy peak in its spectral energy distribution (SED). 
The majority of the blazars detected at very high energies (VHE, E > 100 GeV) by ground- 
based imaging atmospheric Cherenkov telescopes (IACTs) are high-frequency peaked BL 
Lacs (HBL), characterized by synchrotron peaks in the X-ray band (often at energies of 
~ 100 keV). Due to the improved sensitivity of current-generation IACTs such as VERITAS, 
IBLs are attractive targets of observations at VHE gamma rays, particularly because they 
offer the possibility of extension of the VHE blazar catalog to include non-HBLs. VHE 
observations of different blazar classes, including flat-spectrum radio quasars (FSRQs) and 
BL Lac objects, will help in our understanding of the relationship of the different blazar 
populations and, ultimately, the mechanism for particle acceleration and emission in the 
highly-relativistic jets. 



W Com was the first IBL to be detected at very high energies (jAcciari et aj|2008bi ). It 
was discovered as a VHE source by VERITAS during observations carried out over a four 
month period in 2008 (Jan to Apr). During this time a strong gamma-ray outburst was 
measured over a 4-day interval, when the source flared in the middle of March. VERITAS 
reported a steep photon spectrum (r = 3.81 ± 0.35 sta t ± O.34 sys )0 and an integral flux of 
9% of the Crab Nebula flux during the flare nights. The VERITAS detection triggered 
Swift observations, and the multiwavelength data obta ined were adequatel y explained with 



a synchrotron-external Compton (EC) leptonic model (jAcciari et alll2008bl ). 



In this article we report on a second VHE flare in W Com observed by VERITAS in 
2008 June. During this flare, when the source was approximately three times brighter than 
during the 2008 March observations, a multiwavelength campaign was triggered, including 
observations with the space-based AGILE gamma-ray telescope and the Swift and XMM- 
Newton X-ray telescopes. Here we describe the multiwavelength data and derive the spectral 
energy distribution (SED) of the source from contemporaneous data taken throughout the 
flare. 



1 The subscripts stat and sys denote the statistical and systematic error. 
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2. Observations and Results 

A summary of the complete multiwavelength data set on W Com for observation times 
close to the VHE detection on 2008 June 7-9 can be found in Table [1] and Figure [TJ 



2.1. VERITAS: VHE Gamma-ray observations 

VERITAS is an array of four imaging Cherenkov telescopes located at the Fred Lawrence 
Whipple Observatory in southern Arizona. It combines a large effective area (up to 10 5 m 2 ) 
over a wide energy range (100 GeV to 30 TeV) with good energy resolution (15-20%) and 
angular resolution (« 0.1°). The field of view of the VERITAS telescopes is 3.5°. The high 
sensitivity of VERITAS allows the detection of sources with a flux of 1% of the Crab Nebula 
in less than 50 hours of observations. For more details on the VERITAS instrument, see 



Holder et al 


(2006) or 


Weekes et al 


(2002) 



VERITAS observed W Com for 230 minutes on 2008 June 7-9. All observations pass 
quality-selection criteria, which remove data taken during bad weather or affected by hardware- 
related problems. The data were taken in wobble mode, wherein the source was positioned 
at a fixed offset of 0.5° in one of four directions (North, South, East , West) from the c amera 
center. This allows the simultaneous estimate of the backgro und (IFomin et all Il994f) . The 
regions around the VHE gamma-ray blazar 1ES 1218+304 jAcciari et alll2008ch . located 
about 2° north of W Com, and around bright stars (B-band magnitude brighter than 6) are 
excluded from the background estimation. All observations were undertaken in moonlight 
conditions, where the elevated background light levels lead to a lower sensitivity for the 
detection of gamm a rays at the threshold. The threshold of the first-level trigger system 
( iHolder et alll2006l ) was increased to 70 mV (compared to a default value of 50 mV) to allow 
for very high background moonlight levels during observations on 2008 June 9. Table [2] lists 
observation times, elevation range, and background light conditions for the VERITAS ob- 
servations. The different elevations of observation combined with the continuously changing 
background light conditions result in a range of energy thresholds from 200 to 420 GeVo 



The analysis steps consist of calibration and integration of the fla sh- ADC trac es, image 
cleaning, second- moment parameterization of the telescope images (jHillas I Il985l ). stereo- 
scopic re construction of the e vent i mpact position and direction, gamma-hadron separation 



^see e.g. 



Krawczynski et al.l (120061 )). and the generation of photon maps. Most of the far 



2 The energy threshold is defined as the energy at which the peak of the differential counting rate for a 
Crab Nebula-like spectrum occurs. 
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more numerous background events are rejected by comparing the parameterized shape of the 
event images in each telescope with the expected shapes of gamma-ray showers modeled by 
Monte Carlo sim ulations. Me a n-reduc ed- scaled width and mean-reduced- scaled length cuts 
(see definition in lAcciari et all (l2008al )). and an additional cut on the arrival direction of 
the incoming gamma ray (6 2 , defined as the square of the angular distance to the position 
of W Com to the reconstructed shower direction), reject more than 99.9% of successfully 
reconstructed cosmic-ray background events while keeping 45% of the gamma rays. The 
cuts applied here are: integrated charge per image >75 photoelectrons, mean-reduced-scaled 
width and length between -1.2 and 0.5, and G 2 < 0.015 deg 2 . The number of background 
events in the source region are estimated f rom the same field of view using the "reflected- 
region" model with 10 background regions (lAharonian et all 1200 ll ). 

The energy of each event is estimated from detailed Monte Carlo simulations of exten- 
sive air showers and the response of the telescopes, focal plane detectors and electronics. 
The energy reconstruction algorithm uses lookup tables and determines the energy of an 
event as a function of impact parameter, integrated charge per image, background light 
level, offset of the arrival direction from the center of the camera, and zenith and azimuth 
angle. Gamma-ray collection areas for t hese different obsery ing conditions are calculated 
using the same Monte Carlo simulations (IMohanty et allll998l ). The finite energy resolution 
is taken into consideration by calculating collection areas as a function of reconstructed en- 
ergy. The dependence of the collection area on the spectral index is taken into account by 
an iterative process, where collection areas are calculated using the spectral index obtained 
in the previous step. Convergence is usually achieved after 2-3 steps. The spectral recon- 
struction algorithm assigns to each event a collection area according to its estimated energy, 
background light level, offset of the arrival direction from the center of the camera, assumed 
spectral index, and zenith and azimuth angles. Varying conditions, like changing elevations 
or background light levels are therefore taken into account in the flux calculations and spec- 
tral energy reconstruction. It should be noted that the definition of energy threshold used 
here implies that collection areas are non-zero below the stated threshold value. Gamma 
rays are collected, although with lower efficiency, at energies well below 200 GeV even for the 
brightest background light levels. The systematic error in the estimation of the gamma-ray 
energy is dominated by variabilities and uncertainties in the atmospheric conditions, overall 
Cherenkov photon collection efficiency, and limitations of the Monte Carlo simulations. 

Figure [2] shows the sky around W Com as seen by VERITAS in VHE gamma rays. 
A significant flux of very-high-energy gamma rays from W Com is detected by VERITAS 
for the entire data set taken on 2008 June 7-9. A total of 117 excess events (195 on-source 
events and 78 normalized off-source events, normalization factor of 0.10) are measured. This 
corresponds to a significance of 10.3 standard deviations, calculated following Equation 17 
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in lLi &: Mai (119831 ) . Table [2] lists the daily significances and fluxes above 200 GeV, assuming 
a power-law like spectral shape with a photon index of 3.68 (see next paragraph); Figured] 
shows the light curve for these observations. W Com is not detected on 2008 June 9 (MJD 
54626), but observations were restricted to only 32 min due to very high background light 
levels caused by the Moon. The average flux on 2008 J une 7-8 is 2.5-3 tim es higher than 
during the gamma-ray flare from W Com in March 2008 (jAcciari et alll2008bl ). The position 
of the peak of the gamma-ray excess, reconstructed by fitting a 2D-Gaussian function to the 
uncorrelated e xcess sky map, is in agreement with the position of the radio source associated 
with W Com jFev et allhotM ): A RA = 40" ± 31' s ' tat , A dec = -55" ± 41.4^. The systematic 
uncertainty on the pointing, verified with optical pointing monitors, is less than 50". The 
morphology of the excess is compatible with the distribution expected from a point source. 



The differential photon spectrum between 180 GeV and 3 TeV for the measurements 
from 2008 June 7-8 is shown in Figure [3j The shape of the spectrum is consistent with a 
power law dN/dE = C x (E/400 GeV)~ r with an photon index T = 3.68 ± 0.22 stat ± 0.3 
and a flux normalization constant C = (6.5 ± 0.9 s t a t ±1-3 



sys, 



x 10" 11 cm-V^eV- 1 . 



sys 

For 



comparison, the flare in VHE gamma rays from W Com in 2008 March (jAcciari et alll2008bl ) 
is well fit by a power law with a consistent T 
constant C = (2.00 ± 0.31 stat ± 0.5 



3.81 ± 0.35 stat ± 0.34 s „ s , but smaller flux 



sysi 



J sys) 



x 10" 11 cm-V^eV" 



2.2. AGILE: HE Gamma-ray observations 



The Gamma- Ray Imaging Detector (GRID, 30 MeV - 30 GeV) onboard the high energy 
astrophysics satellite AGILE (ITavani et all l2008al) pointed towards W Com continuously 



from 2008 June 9 (18:00 UT) to 15 (12:00 UT) Jverrecchia et all 



2008). The GRID data 



is analyzed using the AGILE standard pipeline ( Vercellone et al I2OO8I ). with a bin size of 
0.25° x 0.25°. Only events flagged as gamma rays and not recorded while the satellite 
crossed the South Atlantic Anomaly are accepted. Events with reconstructed direction less 
than 10° of the Earth limb are rejected, thus reducing contamination from Earth's gamma- 
ray albedo. W Com was observed about 3 degrees off-axis with respect to the boresight and 
a 3.7a excess (pre-trials) of events >100 MeV is found from 12 (03:00 UT) to 13 (03:00) 
June 2008, corresponding to a flux of (90 ± 32) x 10 -8 ph s _1 cm -2 . It should be mentioned 
that this flux is roughly a factor of 1.5 higher t han the highest flux d etected by the Energetic 
Gamma Ray Experiment Telescope (EGRET; lHartman et all (1l999h ) onboard the Compton 
Gamma Ray Observatory and significantly higher than the weekly averaged peak flux of 
(17.2 ± 3.5) x 10~ 8 ph s _1 cm -2 reported by the Large Area Telescope onboard the Ferm i 



Gamma- Ray Space Telescope during its first three months of operation (lAbdo et alll2009bl ). 
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No excess > 3cr is found in the rest of the observing period and upper limits are obtained; 
results can be found in Table [3] and Figure [TJ 



SuperAGILE, the hard X-ray imager onboard AGILE (18-60 keV, iFeroci et all (120071 )) 
observed the source for a net exposure time of 253 ks. The source position in the orthogonal 
SuperAGILE reference system is ~ (3,0) deg, which means that the exposed area is close 
to the full on-axis effective area ( iFeroci et all 120071 ). W Com is not been detected with 
SuperAGILE, and we estimate a 3a upper limit in the 20-60 keV energy of 6 mCrab ~ 
6.9 x 10~ n erg cm~ 2 s _1 (assuming a photon index of F = 2.1). 



2.3. Swift and XMM- Newton: X-ray observations 



Observations of W Com with the Swift satell ite dGehrels et alll2004f) were taken on 2008 
June 7-9. All Swift X-ray telescope (XRT) data (IBurrows et alll2005l ) are reduced using the 
HEAsoft 6.5 package. Event files are calibrated and cleaned following the standard filtering 
criteria using the xrtpipeline task and applying the most recent Swift XRT calibration files. 
All data were taken in Photon Counting (PC) mode, with grades 0-12 selected over the 
energy range 0.3-10 keV. Due to photon pile-up in the core of the point spread function 
(PSF) at rates larger than 0.5 counts s _1 (PC mode), the source events are extracted from 
an annular region with an inner radius of 3 pixels and an outer radius of 30 pixels (47.2 
arcsec). Background counts are extracted from a 40 pixel radius circle in a source-free 
region. Ancillary response files are generated using the xrtmkarf task, with corrections 
applied for the PSF losses and CCD defects. The response matrix Version 11 from the XRT 
calibration files is applied. To ensure valid x 2 minimization statistics during spectral fitting, 
the extracted XRT energy spectra are re-binned to contain a minimum of 20 counts in each 
bin. The spectra can be described by a single power-law convolved with galactic and local 
absorption. Table H] summarizes the observations along with the best fit model parameters. 



W Com was observed by the XMM-Newton Observatory (IJansen et al. Il200ll ) between 
2008 June 14 and June 18 over three consecutive orbits. The three observations comprise data 
from the EPIC detector (0.2-10 ke V) in Small Window mode. The data have been analysed 
using SASv7.1 (IGabriel et all 120041 ) . Several filtering criteria have been applied to the EPIC 
data, including filtering for time periods of high background activity following the standard 
procedure, and filtering only for single- and double-pattern events for EPIC-pn and single 
to quadruple for EPIC-MOS, as well as including only events with good quality (quality 
FLAG=0). For the spectral analysis, circular source and annular background extraction 
regions centered on the source are selected by maximizing the signal-to-noise ratio. The 
spectra are re-binned in order not to oversample the intrinsic energy resolution of the EPIC 
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cameras by a factor not more than 3, while making sure that each spectral channel contains at 
least 25 background-subtracted counts. This allows the use of the \ 2 quality-of-fit estimator 
to find the best fit model. Fits are performed in the 0.2-10 keV energy range simultaneously 
for the three EPIC cameras, where the systematic difference between the EPIC cameras 
is belo w ~5% in nor malization. For the spectral analysis and fitting procedure XSPEC 
vl2.4 (lArnaud Ill996l ) is used. The data can be best described similar to the XRT data by 
a single power-law convolved with galactic and local absorption. Table [5] summarizes the 
observations along with the best fit model parameters. 

The measurements reveal strong variability in X-rays on time scales of much less than 
one day. Fig [1] (panel C) and H] show that the X-ray flux changed by a factor of two during 
the VHE high state on MJD 54625. This i s comparable to observations of W Com with 
BeppoSax in 1998 by iTagliaferri et all (120001 ). where flux variations of a factor of three in 
less than 5 hours is reported. The X-ray flux during the VHE low state of June 2008 is very 
similar to the X-ray activity measured during the detection of W Com in March 2008 (see 
Figure Hj) . 



2.4. Optical, Near-IR, UV and Radio Observations 



Eight optical, one near-IR and three radio observatories contributed data sets to this 
campaign; see Table [1] for an overview. The majority of t he ob servatories are part of the 
GLAST-AGILE Support Program (GASP, see lvillata et~al (boosh ). a subgroup of the Whole 
Earth Blazar Telescope (WEBT§|). In the period considered here, optical observations of W 
Com were carried out at the following observatories: Abastumani, Crimean, Roque de los 
Muchachos (KVA), Talmassons, Torino (for details concerning these observatories, see ref- 
erences provided by WEBT), San Pedro Martir, Northern Optical Telescope (NOlQ) and 
Sapienza University (Italy). Magnitude cal ibration is obtained with respect to the photo- 
metric sequence by iFiorucci fc Tostil (119961 ) . Near-infrared (JHK) data were acquired at 
the AZT-24 tele scope at Campo Im peratore Observatory (Italy). Swift UV/Optical Tele- 
scope (UVOT) (jRoming et all 120051 ) observations were taken in the photometric bands of 
UVW1 (centered at 2600 A), UVM2 (centered at 2246 A), and UVW2 (centered at 1928 A) 
(IPoole et alll2008l ). The uvotsource tool is used to extract counts from the UVOT, correct 
for coincidence losses, apply background subtraction, and calculate the source flux. The 
standard 5 arcsec radius source aperture is used, with a 20 arcsec background region. 



3 http : / / www/oato . inaf . it /blazars / webt 
4 http:/ /www. not. iac.es/ 
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At radio fre quencies, data at 43 GHz were taken with the 32 m antenna at Not e- 
faach et alll2007h . at 14.5 GHz with the 26 m telescope of t he UMRAO jAller et all l20Q3h . 
and at 36.8 GHz with the 13.7 m Metsahovi radio telescope (ITerasranta et allll998l ). 



Data reduction of the optical and radio data followed standard methods and procedures, 
and we refer to the above papers for details. The near infrar ed, optical and UV data are 
corrected for absor ption in our Galaxy using the dust maps of ISchlegel et all (119981 ) and the 
extinction curve of ICardelli et all (119891 ) . Since the blazar is observed in a bright state (see 
Section [3D, a host galaxy contribution has not been subtracted. 



3. Modeling and Discussions 



The single-epoch SEDs for two different time intervals are shown in Fig. [5j The broad- 
band SEDs of W Com show double-humped structures, as found in all known gamma-ray 
blazars. The photon-production mecha n ism in these objects are su ccessfully modeled by 



lepto nic fe.g. IBottcher fc Chiand (12 002) 



Ghisellini fc Madau I (11996^ and hadronic model 



Bottcher et all J2002h . iMiicke et aj ( 2003h . lAharonian I ( 2000p )~The data presented here 



(e.g. 

have been modeled using a leptonic one-zone jet model. F or thi s purpose, a quasi-equilibrium 
version of the model described in IBottcher &: Chiang! (120021 ) is adopted. In this model, 
the observed electromagnetic radiation is interpreted as originating from ultrarelativistic 
electrons (and positrons) in a spherical emission region of co-moving radius Rb, which is 
moving with a relativistic speed /?rc, corresponding to the bulk Lorentz factor T. Lacking 
more detailed constraints on the viewing angle 9 between the jet direction and the line of 
sight, we fix 9 to be the superluminal angle, for which the bulk Lorentz factor V equals 
the Doppler factor D = (T[l — fir cos 9})" 1 , which determines the Doppler shift of photon 
energies and relativistic boosting of intensities. We note that our results mainly depend on 
D so that alternative combinations of V and 9 yielding the same Doppler factor as the ones 
used in our model calculations are also possible, although minor differences in the flux of the 
external- Compton emission with respect to other radiation components (see below) would 
result (see, e.g., Dermer 1995). 

In our calculations, the size of the emission region is constrained by the shortest observed 
variability time scale St VSLTiain through R B < c5t va _ Tjmin D/{1 + z). In the optical and X-rays 
(IBottcher et all 120021 ). variability down to time scales of a few hours has been observed, 
limiting the blob radius to Rb < 10 15 (<5t vai ., m in/hr) (D/10) cm. 



Ultrarelativistic electrons are assumed to be instantaneously accelerated into a power- 
law distribution in electron energy, E e = r ym e c 2 , as Q(j) = Qol~ 9 with a low- and high- 
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energy cutoff 71 and 72, respectively. An equilibrium between this particle injection, radiative 
cooling, and escape of particles from the emission region yields a temporary quasi-equilibrium 
state described by a broken power-law. The particle escape is parameterized through an 
escape time scale parameter 77 > 1 as t esc = rjR/c. The balance between escape and radiative 
cooling will lead to a break in the equilibrium particle distribution at a break Lorentz factor 
7b, where t csc = ^cooi(t)- The cooling time scale is evaluated self-consistently taking into 
account synchrotron, synchrotron-self- Compton (SSC) and external Compton (EC) cooling. 
Depending on whether j b is greater than or less than 7 1; the system will be in the slow cooling 
or fast cooling regime. In the fast cooling regime (7^ < 71), the equilibrium distribution will 
be a broken power-law with 71(7) oc 7~ 2 for 7& < 7 < 71 and 71(7) oc 7~(' 3+1 ) for 71 < 7 < 72- 
In the slow cooling regime (7^ > 71), the equilibrium distribution will be 71(7) oc 7 -9 for 
71 < 7 < 7b and 71(7) oc 7~( ?+1 ) for 75 < 7 < 7 2 . The number density of injected particles is 
normalized to the resulting power in ultrarelativistic electrons propagating along the jet, 



The magnetic field B in the emission region is a free parameter. The corresponding 
Poynting flux along the jet is Lb = 7ri? 2 T 2 /3r cub, with the magnetic energy density ub = 
B 2 /(8n). For each model calculation, the resulting equipartition parameter, e B = L B /L e , is 
evaluated. Modeling results of a large number of blazars, in particular flat-spectrum radio 
quasars, have shown that leptonic models can achieve reasonable fits with the emission region 
being close to equipartition, typically 0.1<,es<,l. However, there is no a priori argument 
which would dictate quasi-equipartition. Therefore, while we disfavor possible fit results 
with e B far from unity, we can not strictly rule out such scenarios. 

Once the quasi-equilibrium particle distribution in the emission region is calculated, our 
code evaluates the radiative output from synchrotron emission, SSC, and EC emission self- 
consistently with the radiative cooling rates. For the EC component, we assume an external 
radiation field which is isotropic in the stationary AGN rest frame and can be approximated 
by a thermal blackbody with peak frequency z/ cxt and radiation energy density u cxt . The 
latter two quantities are free model parameters. The direct emission from this external 
radiation field is added to the emission from the jet to yield the total model SED which we 
fit to the observations. 

In all model calculations the luminosity distance to W Com has been calculated using 
standard A-CDM cosmology with Q m = 0.3 and = 0.7. Absorption of high-energy 
gamma rays by the extragalactic background light is taken into account using the model of 



OC' 




(1) 
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Franceschini et all (120081 ) 



We fit the VERITAS flare detection and high X-ray state (MJD 54624.0 - 54626.0) with 
a pure SSC model, i.e., without any external radiation fields, and with a model with an EC 
component. A Doppler factor of 20 (i.e., D = F = 20) consistent with all observational 
constraints, and well in the range of Doppler factors commonly adopted in other blazar 
modeling works, allowed acceptable fits to the SEDs. We therefore fixed D = F = 20 for the 
remainder of the fitting procedure. 

For a pure SSC fit, the free parameters were thus (1) L e , the injection power of ul- 
trarelativistic electrons into the emission region, (2,3) 71 and 72, the cutoffs of the injected 
electron distribution, (4) q, the injection spectral index, (5) B, the magnetic field, (6), Rb, 
the radius of the emission region, and (7) 77, the particle escape time scale parameter. The 
injection spectral index is tightly constrained by the observed X-ray energy spectral index 
a = q/2, since electrons emitting synchrotron radiation in the X-ray regime are always above 
the critical Lorentz factor 75. The radius of the emission region is constrained through the 
minimum variability time scale of a few hours, as mentioned above. Together with the value 
of the magnetic field, the low-energy cutoff 71, determines the location of the synchrotron 
and gamma-ray peaks in the SED, while the high-energy cutoff 72 influences the location of 
the high-energy cutoffs of the SED, in particular the synchrotron component. The cutoff of 
the SSC component is, in addition, strongly influenced by Klein-Nishina effects. Parameters 
of the SSC fit shown in Fig. [5] are listed in Table El 

No SSC model fit was possible with the emission region being close to equipartition. 
Since there is virtually no observational constraint on the high-energy emission in the low 
(MJD 54626) and intermediate X-ray state (MJD 54631), we could choose a low injection 
power and relatively high magnetic field to achieve a synchrotron peak flux comparable to 
the flaring state, but at a much lower SSC flux. Such a choice of parameters allowed us 
to bring the system close to equipartition. However, almost any positive detection either 
in the Fermi or the VHE gamma-ray range could rule out this interpretation. In the SSC 
interpretation, the most significant difference between the various states consists of a change 
in the electron injection spectral index q from 2.55 in the flaring state to 3.50 and 3.40 in 
the low and intermediate state. 

For a model with an EC component, two more parameters need to be specified: (8) the 
peak frequency i/ ext and the energy density u ext of the external radiation field. As with the 
SSC model, the electron spectral index q is tightly constrained through the X-ray spectral 
index, while the variability time scale constrains the radius of the emission region. In order to 
avoid the problem of required large injection powers (to obtain a high SSC flux) and accord- 
ingly small magnetic fields (not to overpredict the synchrotron flux), the VHE gamma-ray 
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emission can be interpreted as external-Compton emission. In order for Comptonization of 
an external radiation field to be efficient out to gamma-ray energies of E>,Eybe = 300 GeV, 
the external radiation field has to peak at energies E ext <(m e c 2 ) 2 /Eyn E ~ 0.9 eV, i.e., in the 
near-infrared. Therefore, line emission from a putative broad-line region (for which there 
is no evidence in W Comae), would have a too high photon energy characteristic to serve 
efficiently as a source photon field for EC scattering to produce an IC spectrum with peak 
energy near the VHE gamma-ray band. It is therefore more likely that infrared emission, 
e.g., from a near- nuclear dust torus, dominates the external radiation field responsible for 
EC emission at VHE gamma rays. We find that an external radiation field peaking at 
^cxt = 1-5 x 10 14 Hz can, at the same time, serves as an efficient source for EC emission 
and explains the slight near-IR bump in the SED of W Comae. This bump could also be 
due to the host galaxy, and future observations of variability of the IR component or very 
high-resolution imaging are required to break this degeneracy. With the assumption of such 
an external radiation field, acceptable fits to each of the states of W Comae can be achieved 
within a factor of ~ 3 of equipartition. The parameters of our SSC+EC fit are listed in 
Table E 



4. Conclusions 



W Com belongs to the IBL class of blazars, a group with a now-growing number of 
VHE-detected blazars. Other blaz ars detected with V HE gamma rays th at are not of th e 
HBL class include the IBLs 3C 66A (Ucciari et alll2009h and PKS 142 4+240 done et aJl2009h . 
low-frequency-peaked BL Lac objects ( LBLs) such as BL Lacertae ([Albert et all 120071 ). and 
the flat-spectrum radio quasar 3C279 ([Albert et all 120081 ). In this article, we described a 
second VHE flare measured from W Com by VERITAS. The object was detected by VERI- 
TAS at a significance level of 10.3 standard deviations during 2008 June 7-8. The VERITAS 
observations triggered a multiwavelength campaign including AGILE gamma-ray, Swift and 
XMM-Newton X-ray, UV, optical and radio observations. We have carried out extensive 
modeling of the SED of W Com constructed from this contemporaneous multiwavelength 
data set, using a leptonic model considering synchrotron, SSC, and external-Compton emis- 
sion. The SED can be modeled by a simple leptonic SSC model, but the wide separation 
of the peaks in the SED requires a rather low ratio of the magnetic field to electron energy 
density of e# = 2.3 x 10~ 3 . The SSC+EC model returns magnetic field parameters closer 
to equipartition, providing a satisfactory description of the broadband SED. T hese findings 
are sim ilar to the results obtained from the first W Com VHE flare reported by lAcciari et al 
d2008bh . 
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The strong variability of W Com at X-ray and gamma-ray energies on time scales of days 
or less shows that only truly contemporaneous data can provide serious constraints on the 
various emission models. Future observations with VERITAS and the Fermi Gamma-Ray 
Space Telescope should provide even more detailed data to better resolve the short variability 
timescales, helping to further constrain model calculations. 

This research is supported by grants from the US Department of Energy, the US Na- 
tional Science Foundation, and the Smithsonian Institution, by NSERC in Canada, by Sci- 
ence Foundation Ireland, and by STFC in the UK. We acknowledge the excellent work of the 
technical support staff at the FLWO and the collaborating institutions in the construction 
and operation of the instrument. Financial support by the Italian Space Agency through 
contract ASI-INAF 1/088/06/0 is acknowledged. Support of UMRAO from NSF and Uni- 
versity of Michigan is acknowledged. NVE and VML acknowledge support from RFBR 
grant 09-02-00092. The GASP-WEBT consortium is acknowledged. We acknowledge the 
efforts of the Swift team for providing the UVOT and XRT observations. This work was 
partially supported by NASA through XMM- Newton Guest Observer Program award No. 
NNX08AD67G and the Swift Guest Investigator Program award No. NNX08AU13G. The 
Metsahovi team acknowledges the support from the Academy of Finland. 

Facilities: VERITAS, AGILE, Swift, XMM-Newton, GASP-WEBT 

REFERENCES 

Abdo, A.A. et al (Fermi LAT Collaboration) 2009a, ApJS, 183, 46 
Abdo, A.A. et al (Fermi LAT Collaboration) 2009b, ApJ, 700, 597 
Acciari, V.A. et al (VERITAS Collaboration) 2008a, ApJ, 679, 1427 
Acciari, V.A. et al (VERITAS Collaboration) 2008b, ApJ, 684, L73 
Acciari, V.A. et al (VERITAS Collaboration) 2008c, ApJ, 695, 1370 
Acciari, V.A. et al (VERITAS Collaboration) 2009, ApJ, 693, L104 
Aharonian, F. 2000, New Astronomy, 5, 377 

Aharonian, F. et al (HEGRA Collaboration) 2001, A&A, 370, 112 
Albert, J. et al (MAGIC Collaboration) 2007, ApJ, 666, L17 



-16- 



Albert, J. et al (MAGIC Collaboration) 2008, Science, 320, 1752 
Aller, M.F. et al 2003, ApJ, 586, 33 

Arnaud, K. A. 1996; in ASP Conf. Ser. 101: Astronomical Data Analysis Software and 
Systems (ADASS) V, 17. 

Bach, U. et al 2007, A&A, 464, 175 

Berge, D. et al. 2007, A&A, 466, 1219 

Bdttcher, M. & Chiang, J. 2002, ApJ, 581, 127 

Bottcher, M., Mukherjee, R. & A.Reimer 2002a, ApJ, 581, 143 

Burrows D. et al 2005, SSRv, 120, 165 

Cardelli, J.A., Clayton, G.C., & Mathis, J.S. 1989, ApJ, 345, 245 

Dermer, CD. 1995, ApJ, 446 L63 

Dickney, J.M. & Lockman F.J. 1990; ARA&A. 28, 215 

Feroci, M., et al. 2007, Nucl. Instrum. Methods Phys. Res. A, 581, 728 

Fey, A.L. et al. 2004, AJ, 127, 3587 

Fiorucci, M. & Tosti, G. 1996, A&AS, 116, 403 

Fomin, V.P. et al 1999, Astroparticle Physics 2, 137 

Franceschini, A. & Rodighiero, G. & Vaccari, 2008, A&A487, 837 

Gabriel, C, Denby, M., Fyfe, D. J. et al. 2004; in ASP Conf. Ser. 314: Astronomical Data 
Analysis Software and Systems (ADASS) XIII, 759. 

Gehrels, N. et al 2004, ApJ, 611, 1005 

Ghisellini, G. & Madau, P. 1996, MNRAS280, 67 

Hartman, R.C., et al. 1999, ApJS, 123, 79 

Helene, O. 1983, Nuclear Instruments & Methods 212, 319 

Hillas M. 1985, Proc. of the 19th ICRC (La Jolla, USA), 3, 445 

Holder, J. et al (VERITAS Collaboration) 2006, Astroparticle Physics 25, 391 



Jansen, F. et al. 2001, A&A, 365, LI. 

Krawczynski H. et al. 2006, Astroparticle Physics, 25, 380 

Li, T.P. & Ma, Y.Q. 1983, ApJ, 272, 317 

Miicke, A. et al 2003, Astroparticle Physics 18, 593 

Mohanty, G. et al. 1998, Astroparticle Physics, 9, 15 

Nieppola, E., Tornikoski, M. & Valtaoja, E. 2006 A&A, 445, 441 

Ong, R.A. (VERITAS Collaboration) 2009, ATEL 2084 

Poole T.S. et al 2008, MNRAS, 383, 627 

Roming, P.W. et al 2005, Space Science Reviews, 120, 95 

Schlegel, D.J., Finkbeiner, D.P, & Davis, M. 1998, ApJ, 500, 525 

Tagliaferri, G. et al. 2000., A&A, 354, 431 

Tavani, M., et al. 2008, Nucl. Instrum. Methods Phys. Res. A, 588, 52 

Terasranta, H. et al, A&AS, 132, 305 

Vercellone, S., et al. 2008, ApJL, 676, 13 

Verrecchia, F., et al. 2008, The Astronomer's Telegram #1582 

Villata, M. et al. 2008, A&A, 481,L79 

Weekes, T.C. et al 2002, Astroparticle Physics 17, 221 

Wing, R.F. 1973, AJ, 78, 684 



This preprint was prepared with the AAS IATgX macros v5.2. 



-18- 



Table 1. Observatories contributing to the presented data set. 



Waveband 


Observatory 


Frequency/Band/ 
Energy Range 


MJD Range 


Radio 


UMRAO 


14.5 GHz 


54630-54633 




Metsahovi 


36.8 GHz 


54623-54634 




Noto 


43 GHz 


54611 


NIR/Optical/UV 


NOT 


U/B/V/R/I 


54636 




rp 1 

luorla 


K 


E A COO E A d A E 

o4ozz - o4o4o 




Abastumani 


R 


54617 - 54637 




Sapienza University 


R 


54627 - 54634 




San Pedro Martir 


R 


54620 - 54624 




KVA 


R 


54626 - 54633 




Crimean 


R 


54623 - 54627 




Talmassons 


R 


54628 




Torino 


R 


54630 




Campo Imperatore 


J/H/K 


54627 - 54633 




Swift UVOT 


U/B/V/UV 


54625 - 54626 


X-ray 


Swift XRT 


0.3-10 keV 


54625 - 54626 




SuperAGILE 


20-60 keV 


54626 - 54630 




XMM-Newton EPIC 


0.2-10 keV 


54631 - 54635 


HE Gamma-ray 


AGILE GRID 


30 MeV - 30 GeV 


54626 - 54632 


VHE Gamma-ray 


VERITAS 


0.1-30 TeV 


54624 - 54626 
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Table 2. Details of VERITAS observations of W Com on 2008 June 7-9. The energy 
threshold for fluxes and upper flux limits (99% confidence level; assuming a photon index 
of T = 3.68) is 200 GeV. Errors are given at the la level. 



MJD 


elevation 


observation 


average pedestal 


significance 


flux or 




range 


time 


variations 3. 


(pre-trials) 


upper flux limit 






[min] 


[dc] 


M 


[cm -2 s _1 ] 


54624.16 - 54624.23 


53-73° 


100.2 


7.8-8.0 


8.9 


(5.0 ±0.8) x 10- 11 


54625.17 - 54625.24 


49-68° 


100.2 


8.1-9.7 


7.9 


(6.2 ± 1.2) x 10~ n 


54626.18 - 54626.20 


59-60° 


32.0 


12.2-12.3 b 


-1.0 


< 3.21 x 10~ n 



a The average pedestal variation in digital counts [dc] indicates the background light level. Values 
of 6.5 to 6.8 are typical for regular observations of extragalactic targets on moonless nights. All 
observations presented here are taken in moonlight conditions. 

b Data taken with increased pixel (PMT) trigger threshold (at 70 mV CFD trigger threshold instead 
of the regular 50 mV). 



Table 3. Details and results of the AGILE GRID observations of W Com on 2008 June 
9-15. The energy threshold for fluxes and upper flux limits (99% confidence level; assuming 
a photon index of T = 2.1) is 100 MeV. Errors are given at the la level. 



MJD 


significance 


flux or 




(pre-trials) 


upper flux limits 






[cm- 2 s- 1 ] 


54626.75 - 54629.12 


< 3(7 


< 60 x 10~ 8 


54629.12 - 54630.12 


3.7a 


(90 ± 34) x 10~ 8 


54630.12 - 54632.50 


< 3a 


< 55.5 x 10~ 8 
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Table 4. Details and results of the Swift/XRT observations of W Com 2008 June 7-9. 
The galactic N# Ga ; has been fixed to a value of 1.88-10 20 cm -2 . The redshift of the source 



was assumed to be 0.102. 


Errors are given 


at the la level. 


MJD 


Exposure 


Photon index 


Flux F 2 _iokeV 




[ksec] 


r 


[10~ 12 ergs cm -2 s _1 ] 


54624.97 - 54624.98 


0.52 


2.49 ±0.19 


3.90 ±0.97 


54625.04 - 54625.05 


0.84 


2.71 ±0.15 


3.70 ±0.76 


54625.11 - 54625.12 


1.38 


2.55 ±0.09 


4.75 ±0.55 


54625.17 - 54625.20 


2.51 


2.36 ±0.05 


9.33 ±0.74 


54625.24 - 54625.27 


2.47 


2.59 ±0.07 


4.62 ±0.37 


54626.11 - 54626.21 


5.07 


2.69 ±0.10 


1.00 ±0.18 



Table 5. Details and results of the XMM-Newton observations of W Com 2008 June 
14-18. The galactic Nhv^ has been fixed to a value of 1.88-10 20 cm -2 as obtained from 



Dickney fc Lockman I (119901 ). The redshift of the source was assumed to be 0.102. Errors 





are g 


iven at the 1 


sigma level. 




MJD 


Exposure 
[ksec] 


[10 20 cm" 2 ] 


Photon index 

r 


Flux F 2 _l0keV 

[10 -12 ergs cm -2 s" 1 ] 


54631.50 - 54631.55 
54633.15 - 54633.17 
54635.14 - 54635.16 


28.0 
16.0 
11.0 


9 9 n +0.09 
z - zu -0.09 
1 QO+0.14 
± - oy -0.13 
i ntr+0.16 
i - ud -0.15 


o 70+0.01 
z - ' a -0.01 

9 00+0.02 
^•oS-0.02 
9 77+0.02 
z - ' ' -0.02 


2 69+ - 02 

z - U3 -0.02 
1 t-q+0.03 

1 89+ - 03 
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Table 6: Parameters of SSC and SSC+EC fits to the SEDs of W Com on MJD 54624.0 - 
54626.0. 



Parameter 


Symbol 


SSC 


SSC+EC 


Doppler factor 


D 


20 


20 


Electron power [erg s~ 4 ] 


L e 


3.4 x 10 44 


5.7 x 10 43 


Blob radius [cm] 


Rb 


3 x 10 15 


1Q 16 


Low-energy cutoff 


7i 


9 x 10 3 


8 x 10 3 


High-energy cutoff 


72 


2.5 x 10 5 


3 x 10 5 


Electron injection index 


Q 


2.55 


2.55 


Magnetic field [G] 


B 


0.24 


0.35 


B-field equipartition parameter 


eB 


2.3 x 10~ 3 


0.32 


Electron escape time scale parameter 


V 


300 


300 


Minimum variability time scale [hr] 


^var,min 


1.5 


5.1 


External radiation peak frequency [Hz] 


^ext 




1.5 x 10 14 


External radiation energy density [erg cm -3 ] 


^ext 




2.4 x 10" 4 
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Fig. 1.— Multiwavelength light curve of W Com for MJD 54622 to 54636. Panel a: VHE 
gamma-ray light curve (E>200 GeV) as measured by VERITAS. The flux in VHE gamma 
rays corresponds to approximately 25% of the flux of the Crab Nebula above 200 GeV. 
Panel b: Gamma-ray light curve (E> 100 MeV) as measured by AGILE. Panel c: X-ray 
(Swift XRT: 2-10 keV; circles; XMM-Newton EPIC: squares) Panel d: Swift UVOT (UVW1: 
squares; UVM2: downward-pointing triangles; UVW2: upward-pointing triangles) Panel e: 
Light curves of negative optical magnitudes (R-Band; filled circles: Tuorla; filled squares: 
Abastumani; filled triangles: San Pedro Martir; diamonds: Sapienza University; open circles: 
KVA; open squares: Crimean; open stars: NOT; open triangles: Torino; open crosses: 
Talmassons ). Panel f: Radio light curve (circles: UMRAO 14.5 GHz; triangles: M etsahovi 



37 G Hz). Downward pointing arrows indicate upper flux limits (99% confidence level; iHelene 
jl983h ). 
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Fig. 2. — Sky map of significances of gamma-ray emission from the region around W Com. 
The background is estimated using the reflected region model (10 backgrou nd regions , over- 
sampling radius 0.12°). The position of W Com derived from radio data ( IFey et all 120041 ) 
is indicated by a white cross. The dashed circles indicate positions of bright stars and 
their B-band magnitudes in the field of view; regions around these stars are excluded from 
the background es timation. Two sources listed in the Fermi bright gamma-ray source list 
([Abdo et alll2009af) . and firmly associated with the blazars W Com and B2 1215 (I Abdo et al 
2009bh . are shown with their 95% confidence £1X6 cl cLS circles with V in their centre. The 



circle at the bottom right indicates the angular resolution of the VERITAS observations. 
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Fig. 3.— Differential VHE photon spectrum for W Com for MJD 54624.16 to 54625.24 (2008 
June 7-8). The markers indicate measured data points and the continuous line a fit assuming 
a power-law distribution. Error bars show statistical errors only. For compariso n, the photon 
spectr um of W Com derived from VERITAS measurements in March 2008 (lAcciari et al 
2008bl ) is indicated by a dashed line. 
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Fig. 4. — Spectral energy distribution for five different X-ray measurements with Swift XRT 
and XMM-Newton EPIC for 2008 June 7-18 and the Swift XRT measurements in March 
2008 (MJD 54539.4). 
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Fig. 5. — High-state spectral energy distribution of W Com for MJD 54624 to 54626 including 
VERITAS, Swift XRT/UVOT, optical and radio data (filled circular markers). The averages 
of the optical, NIR and radio fluxes calculated over the time range from MJD 54610 to 54645 
are shown here . Downward-pointing arrows indicate upper flux limits (99% confidence level) 
faelene Ifl9~83h . For comparison, the VERITAS, AGILE and Swift XRT/UVOT data for 
MJD 54626 to 54626.9 are shown as grey open squares and grey downward pointing arrows. 
Results from synchrotron-self-Compton (SSC) and external-Compton (SSC+EC) models 
are shown as continuous and dashed lines, respectively (see text for details). The different 
components (disk, sync=synchrotron, ec, ssc) of the SSC+EC models are indicated by dotted 
lines. 



